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Abstract. Effects of the surface exchange anisotropy on ordering of ferro- 
magnetic films are studied for the exactly solvable classical spin-vector model 
with D — > oo components. For small surface anisotropy t)j < 1 (defined rela- 
tive to the exchange interaction), the shift of T c in a film consisting of W > 1 
layers behaves as T c bulk - T C (N) oc (1/iV) hx(l/rj' a ) in three dimensions. The 
finite-size-scaling limit T c bulk - T C (N) oc l/(ri' 1/2 N 2 ), which is realized for the 
model with a bulk anisotropy rj <C 1 in the range Nr/' 1 ^ 2 J> 1, never appears 
for the model with the pure surface anisotropy. Here for N exp( — l/r]' s ) !> 1 in 
three dimensions, film orders at a temperature above T^ ulk (the surface phase 
transition). In the semi-infinite geometry, the surface phase transition occurs 
for whatever small values of r]' B (i.e., the special phase transition corresponds 
to Tc ) in dimensions three and lower. 



1. Introduction 

Reduction of Curie temperatures T c of ferromagnetic films consisting of A ^ 1 layers 
with respect to the bulk value is usually represented in the form 

pvbulk _ T c (A0]/T c bulk ~ A/N x . (1.1) 

For the exponent A the finite-size scaling theory |l], || yields A = l/vt, where uj, is 
critical index for the bulk correlation length. The above result has been derived with 
the Ising model, or the field model with one-component order parameter, in mind. 
For weakly anisotropic Heisenberg model and, in general, for models with several 
spin components, the nearly Goldstone modes can drastically change the character 
of ordering in magnetic films. In particular, in the dimensionality range d < 3 in the 
isotropic limit at low temperatures, the film behaves as a system of dimensionality 
d! = d— 1 < 2 and cannot or der because of long- wavelength fluctuations. This means 
that the amplitude A in Eq. (]l.l[) should diverge in the isotropic limit. Moreover, even 



the functional form of Eq. (hi) should change to reflect explicitly the d! dimensional 
nature of a nearly isotropic film. For the model with the uniaxial exchange anisotropy 
(longitudinal spin components coupled by J and transverse components coupled by 
77 J with r\ < 1, so that rj' = 1 — r\ measures the anisotropy) it was shown in Rcfs. 
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|^] that Eq. (1.1) is only valid for rather thick films, Nk c > 1, where n c = l/£ c j_ = 
^2(2(1/77-1) is the inverse transverse correlation length at the bulk critical point, 
which goes to zero in the isotropic limit. Here in three dimensions for the classical 
spin- vector model with D — > oo components one has A = 2 and A ~ 1/k c . In the 
range Nk c <1 a rf' dimensional behavior is realized, which is characterized by A = 1 
and A ~ ln[l/(«; c iV)] in three dimensions, d = 3. For extremely small anisotropics, 
the film orders at such low temperatures that spins along the direction perpendicular 
to the surface are strongly correlated with each other and they can be considered as 
single composite spins. Thus the film is mapped on d'-dimensional monolayer with 
the exchange interaction NJ, which yields 0] 

T c (d,J,rj',N)^T c (d',NJ,dr)'/d',l). (1.2) 

Although the results of Refs. || H have been obtained for the infinite-component 
classical vector model, the qualitative features of this solution should be shared by 
the more realistic Heisenberg model, D = 3. In particular, formula ( |l.2| ) is model 
independent and valid for all D > 2. 

The purpose of this paper is to study ordering in magnetic films with a surface 
anisotropy. The latter arises, typically, due to the violation of the symmetry of a 
crystal field acting on the magnetic ions at the surface. Although this anisotropy has 
a single-site form, we will consider here the anisotropy of the exchange interactions 
between the surfaces spins, instead. This leads to the same qualitative results and 
allows one to use the formalism developed for the exchange- anisotropy models in Refs. 
||, ||, [| [|. One can expect that the surface anisotropy stabilizes ordering in films 
with d < 3 and N 3> 1 weaker than the bulk one. If the surface anisotropy is very 
small, then T c <§; T^k, and at such low temperature, its influence should redistribute 
over N layers, so that its effective value be rj' eS ~ T)' s /N. The latter should result in a 
more pronounced suppression of T c in magnetic films. This can be immediately seen 
in the case of extremely small surface anisotropy, where the analogue of Eq. ( |1.2| ) 
reads 

T c (d, J, V ' S ,N) = T c (d', NJ, rfjN, 1). (1.3) 

A specific feature of the model with pure surface anisotropy is the absence of the 
finite length scale, such as the transverse correlation length £ c .j_ = 1/n, at criticality. 
As a result, the system is always in the range Nk <C 1 and there is no crossover to 
the finite-size-scaling regime of Eq. (^^). As a result, for a small surface anisotropy, 
the corresponding analytical solution for the Curie temperature of the film holds in a 
much wider range of N. 

If surface anisotropy exceeds a critical value r] s ^ c (N), the Curie temperature of the 
film exceeds the bulk Curie temperature: T c > T^ ulk . The possibility of this effect, 
which is absent in the mean field approximation (MFA) , can be seen from the following 
simple arguments. The isotropic large- D model orders at T^ ulk = Jo/(DWd), where 
Jo is the zero Fourier component of the exchange interaction and W ( i = -Pd(l) [see 



Eq. ( 2.21 )1 is t ne Watson integral containing the information on the lattice dimen- 
sionality and structure. On the other hand, the Curie temperature of the monolayer 
with the (surface) anisotropy of the extreme Ising type, ~q s = (i.e., rj' s = 1), is 
T c (l) = (d' /d)Jo/D. For the simple cubic lattice one has W3 = 1.51639, so that the 
Curie temperature of the anisotropic monolayer slightly exceeds the isotropic bulk 
Curie temperature. That is, the lack of interacting neighbours at the surface can be 
compensated for by a stronger suppression of T c bulk due to long-wavelength fluctua- 
tions making contribution to Wd- It is clear that the bilayer has a substantially higher 
value of T c than the monolayer, and that in dimensions lower than 3 the bulk Curie 
temperature is suppressed even stronger. For the continuous-dimension model intro- 
duced in Ref. 0, one has W3.0 = 1.719324 and W 2 . 5 = 2.527059. In two dimensions 



2 



and below,, Wd diverges and thus T^ ulk goes to zero in the isotropic limit. On the 
other hand, the theory predicts a finite-temperature surface phase transition for any 
nonzero values of the surface anisotropy rj' s . Thus, for d < 2 the surface anisotropy is 
the only source of ordering. This situation is realized only in the limit D — > oo, how- 
ever. Since in two dimensions the surface is one dimensional, ordering at the surface 
should be destroyed by thermal fluctuations of the longitudinal spin components for 
any finite D. In fact, surface anisotropy plays a major role already for d < 3. We will 
see below that rj StC (N) goes to zero in the limit N — > oo in this dimensionality range. 
Thus in the semi-infinite geometry a surface phase transition above T^ ulk occurs for 
whatever small value of the surface anisotropy, i.e., the bulk Curie temperature is the 
temperature of the special phase transition as well! 

The main part of this paper is organized as follows. In Sec. ||]the closed system of 
equations describing the D — > oo component spin- vector model in the symmetric phase 
is written down. In Sec. |3] the analytical calculation of the correction to T c in films 
with a weak surface anisotropy is presented. In Sec. |4] the surface phase transition is 
considered. The results of numerical calculations are at appropriate places in sections 
above. In Sec. [5] the results obtained are summarized, and possibilities of finding 
similar regimes in more realistic models are discussed. 



2. Basic equations and their solution 

The Hamiltonian of the anisotropic classical /^-component spin-vector model can be 
written in the form 



H 



- ^ ^ Jjj f m zi m Z j + rjjj ^ m ai m a j j , |mj| = 1, (2.1) 



where dimensionless anisotropy factors satisfy rfcj < 1. This model was introduced, in 
the isotropic form, by Stanley, who showed that its partition function in the spatially 
homogeneous case in the limit D — > 00 coincides with that of the spherical model 
M. There are, however, a number of essential differences between the exactly solvable 



limit D — ► 00 of Eq. (2.1) and the spherical model. In particular, there is only 
one correlation function (CF) in the spherical model, and thus this model cannot 
incorporate anisotropy. In the D — * 00 model, there are longitudinal and transverse 
CFs which differ below T c , even in the spatially homogeneous isotropic case ||. 

The system of equations describing the spatially inhomogcncous D — > 00 model 
both above and below T c was obtained in Ref. At or above T c in zero field, 
the magnetization (m^) is zero and the model is described by the closed system 
of equations for the correlation functions of transverse (a > 2) spin components, 
Sij = D(m a im a j), and the spatially varying gap parameter, Gi. (The definition 
of Gi can be found in Ref. [||; here it is nonessential.) In the film geometry, it is 
convenient to use the Fourier representation in d! = d — 1 translationally invariant 
dimensions parallel to the surface and the site representation in the dth dimension. 
The equations can be easily generalized for the anisotropy factors taking the values 
rj nn within the nth layer and rj n ,n±i for the interaction between the nth and (n± l)th 
layers. For the model with nearest-neighbor (nn) interactions, the equation for the 
Fourier-transformed CF a nn > (q) then takes the form of a system of second-order 
finite-difference equations in the set of layers n = 1, 2, N: 

where b n is given by 

b n = d/( Vnn G n ) - d'X' q , (2.3) 
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Aq for the d-dimensional hypercubic lattice reads 

1 d ' 



d> 

i=l 



and the lattice spacing has been set to unity. In Eq. (2.2), 9 is the reduced temperature 
defined by 

-* —iMFA/-\ _ ^0 



where for hypercubic lattices Jo = 2dJ. The quantities <7o,n' and crjv+i,n' m the 
nonexisting layers, which enter equations (2.2) at the film boundaries n = 1 and 
n = N, are set to 

C0,n' = CJV+l.n' = (2.6) 

as free boundary conditions. The autocorrelation functions in each of N layers, s nn , 
satisfy the set of constraint equations 

f d d 'q 

which are the consequence of the spin rigidity |nij| = 1. A straightforward algorithm 
for numerical solving the equations above is to compute, for a given set of G n , all a nn 
from the system of linear equations (2.2) and then insert the results in Eq. (2.7) to 
obtain, after the integration over the Brillouin zone, a set of nonlinear equations for 
G n . 

The first step of the routine described above can be conveniently done with the 
help of the continued-fraction formalism which is described in detail in Refs. || |J. 
For a particular type of the model with surface and bulk anisotropics, which is defined 
by 

7/11 = VNN = Vs < 1, V = Vn,n±l = V < 1 11,NN), (2.8) 



and which will be studied below, it is convenient to rewrite equations (2.2) in the 
form 

26„<7„„' " OVi+l,n' - Vn-l,n' = ( 2c $ /v)$nn' , (2.9) 

where b n — (r] s /ri)b n for n = 1, N and b n — b n otherwise. Explicitly, 

K = d/(r)G n ) - rf'Aq + (1 - r) s /r))d'\'q(6 n>1 + S nN ). (2.10) 

An alternative way to find er lm < , which is more appropriate for the analytical 
treatment, is to represent equations (2.2) in the matrix form 

Ba = diag(2d0/?7 nn ), B nn = 2b n , S n ,T.±i = -»?n,n±i/V) ( 2 - n ) 

so that the solution for a nn ' is given by cr nn / = (2d9/r] n i n i)B~^, . Since the diagonal 
part of the matrix B, which depends on the wave vector q, is proportional to the 
unity matrix, the eigenvalues and eigenvectors of B can be defined as 

BU p = [n P + 2d'(l-\' q )]U p , p= 1,2,..., N, (2.12) 

the eigenvectors U p being independent of q. It should be noted that matrix B is 
nonsymmetric, B„, n ±i = -t)n,n±i/v ^ B n ±x <n = -?7n,n±i/??n±l,n±i, if anisotropy 
factors r)nn change from one layer to the other. In this case its left eigenvectors Wj 
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differ from its right eigenvectors U p . The Green function cr„„< can be expanded over 
the set of eigenvectors of the problem as follows 



2d9 "U^, 



Here matrix U is composed of the right eigenvectors U p as columns and W T is com- 
posed of the left eigenvectors WT as raws. The right and left eigenvectors satisfy the 
biorthogonality condition ^2 n W np U n pi = 5 pp i. In general, matrix U in nonunitary: 
If- 1 = VK T 7^ £> T . Integration in Eq. (§/7|) can be performed analytically with the 
result 

where Pd'(X) is the lattice Green function for the layer, which is defined similarly to 
the lattice Green function P(X) = Pd(X) below. Using this method with tabulated 
values of Pd> (X) can save computer time, in comparison to the continued-fraction 
method. On the other hand, the continued-fraction method is fast enough and already 
implemented, so that it will be used here. The diagonalization formalism above will 
be used for analytically solving the problem in the next two sections. 

After the set of G n for a given temperature has been determined, one can compute 



the longitudinal CF a^,(q) from Eqs. (2.2) and (|2.3|), where all anisotropy factors 



r] nn > are replaced by 1. The Curie temperature of the film 9 C can now be found from 
the equation 

[<C(q = 0)]- x =0. (2.15) 

In a usual situation, the above condition should be used in the middle of the film, 
n ~ N/2, because for large N the critical divergence of the spin CF at the surface 
is suppressed [[l0| [ll], If ordering of the film is driven by the surface, it is more 



convenient to use Eq. ( |2 .15 ) for n = 1. This equation has, in general, N roots, as 
we will see below. One should choose the maximal root for 9 C , all other roots are 
unphysical. Below 9 C , the spontaneous magnetization appears, and the very form of 
the equations change. 



One can also represent cr zz n , in the form of Eq. ( 2.13Q with r) n > n > =>■ 1, where 



eigenvalues \ z p and eigenvectors components U z p correspond to the problem with the 



matrix B z . The latter is defined by Eq. ( 2.11 ), where anisotropy factors T) nn > are 



replaced by 1. Since B z is a real symmetric matrix, {B Z ) T — B z , matrix U z is 
unitary: (U 2 )^ 1 = (U Z ) T , i.e., U z ^ x = U z p . The eigenvalue problem corresponding 
to the longitudinal CF can be written in the form of a discrete Schrodinger equation 
for a particle with mass m = 1/2: 

-iP n - 1 +2iP n -^ n+1 + V n ^ n = Eip n , V n =2d(l/G n -1), (2.16) 

as in quantum tight-binding models. This form is useful for the interpretation of the 



results; the eigenvectors and eigenvalues of Eq. (2.16) are more compact forms of the 
quantities introduced above: 

o£, (q)=2d^|^, <?«!, (2.17) 

where E p = [i z p and ip np = U z p . The condition for the Curie temperature of the film 
has the form 

E 1 (0 C )=Q, (2.18) 
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where E\ is the lowest of the eigenvalues E p . The N — 1 solutions corresponding to 
E p = 0, p > 2, are unphysical. It should be noted that for the transverse correla- 
tion function the problem cannot, in general, be interpreted quantum mechanically, 
since the matrix B may be non-Hermitean, as is the case for the model with surface 
anisotropy. The eigenvalues of the transverse problem, fi p , exceed the longitudinal 
eigenvalues E p ; in the Ising limit r\ <C 1 one has /x p cx 1/f], whereas E p become 
independent of r\. 

One should note that the longitudinal CF is in our formalism only a "slave" 
quantity, it does not affect the basic equations of the model and is not subject to a 



constraint condition similar to Eq. ([2/?]). The physical reason for that is irrelevance 
of fluctuations of the single longitudinal component in comparison to those of D — 1 
transverse ones in the limit D — > oo. 

In the spatially homogeneous bulk sample one has G n = G and rj nn = T] n ,n±i = V: 
and the transverse CF can be easily found || |J] . The resulting equation for the gap 
parameter has the form 

0GP(r]G) = 1, (2.19) 

where 

f d d V 1 

P(X)= - (2.20) 

is the lattice Green function. The quantity A^. = J^/Jo for the nearest-neighbour 



interaction is given by Eq. (2.4) with d! d and q =>• k. The solution G of Eq. 
( 2.19| ) increases with lowering temperature 6; at G = 1 the gap in the longitudinal 



CF closes, longitudinal susceptibility diverges, and the phase transition occurs. This 
defines the bulk transition temperature [fiT 



nbulk 



1/P( V ), (2.21) 



that generalizes the well known result for the spherical model 8 C = l/P(l) S. The 
lattice Green function P(X) satisfies P(0) = 1 and has a singularity at X — > 1, the 
form of which in different dimensions can be found in Ref. |@|. For d < 2, the Watson 
integral W = P(l) goes to infinity; thus formula ( 2.2l| ) yields nonzero values of the 



Curie temperature only for the anisotropic model, r] < 1. It should be noted that 
in the anisotropic case the critical indices of the model coincide with the mean-field 
ones due to the suppression of the singularity of P(rjG) for G — » 1. Below C , the 
spontaneous magnetization appears, and G sticks to 1. 



In Eq. ( 2.2C ) one has = 1 — k /(2c?) in the long- w avele ngth limit. Thus the 



inverse transverse correlation length n following from Eq. ( |2.20 ) is defined by 



k 2 = 2d[\/{rjG) — 1]. (2.22) 



Its critical-point value k c = y/2d[l/rj — 1] measures the bulk anisotropy and varies 
between for the isotropic model and oo for the classical Ising model. The inverse 
longitudinal correlation length n z is determined by n 2 z = 2d[l/G— 1] and it diverges at 
the critical point. In contrast to finite- -D theories, where the longitudinal correlation 
length £ cz = 1/k z plays the major role in the scaling, here in the limit D — > oo it 
becomes only a slave variable, whereas all the physical quantities, except the longitu- 
dinal CF, are scaled with the transverse correlation length £ C) j_ = 1/k ^ |j. 

Numerical solution of the problem with the method described in this section above 
yields the results for 9 C (N) of the three-dimensional film with a simple cubic structure, 
which are shown in Fig. |l|. One can see that for small transverse anisotropies, C (N) 
approaches its bulk limit much slowlyer than the solution for the model with bulk 
anisotropy i] — (classical Ising model), which is shown by solid circles. Since in the 
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Fig. 1: Curie temperatures of the iV-layer simple-cubic- lattice film for different values 
of surface anisotropy. Horizontal dotted line is the bulk value of T c , Solid circles are 
the values of T c for the model with the bulk anisotropy 77' = 1 (classical Ising model) . 



latter case transverse spin components are switched off, the result coincides with that 
of the mean field approximation [ fi"3| 

W) = l-Ul-cos^W-l(£) 2 (N»l). (2.23) 



As can be seen from Fig. [l], the suppression of the Curie temperature in films 
with a weak surface anisotropy may be quite pronounced, especially in comparison 
with the mean-field result shown by solid circles. Large T c shifts in films have been 
observed in many experiments (see, e.g., Ref. ]f4[). For larger values of rf a , the film 
orders above T^ ulk . This is an indication that ordering at the surface occures first and 
thus determines the Curie temperature of the film. The decrease of 9 C (N) with N in 
this region can be easily explained. For N — ► 00 both surfaces order independently 
at some 8 c (oo). For finite N, the surfaces interact with each other across the film 
and thus help each other to order. The interaction between surfaces, and thus the 
corresponding increase of 9 C , should decay exponentially with the film thicknes N, 
the characteristic length being the bulk correlation length. One can see that for 
■q' s = 1 surfaces order at a temperature substantially higher than 8^ ulk , where the 
bulk correlation length is rather short. With lowering rj' s , the bulk correlation length 
at 9 C increases, and the effect of the interaction of surfaces becomes more and more 
pronounced. The mechanism described above will be considered in more detail in Sec. 

I 

The Curie temperatures of films consisting of one and two layers can be calculated 
analytically since there is no inhomogeneity of the gap parameter G n . For the mono- 



layer the result can be obtained by a straightforward renormalization of Eq. (2.21) 
and has the form Q^ 1 = [d/(d — l)]Pd'(r]) (there is no difference between the models 
with bulk and surface anisotropies). For the bilayer, the surface-anisotropy model 
orders, evidently, at lower temperatures than the bulk-anisotropy one. For the latter, 
the expression for 9 C can be found in Ref. Q . For the surface-anisotropy model, the 
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Fig. 2: Curie-temperature shifts in simple-cubic-lattice films for small values of surface 
anisotropy. Dashed lines represent Eq. ( 2.25|) with C3 = 1.36. 



result has the form 



d! 



(2.24) 



For r] s = one has P = 1, and this formula yields 8 C — 2(d — I)/ (2d — 1), which 
becomes 4/5 for d = 3 (see Fig. [|). For com parison, for the model with the bulk 
anisotropy 77 = 1, the mean-field formula ( |2.23 ) yields 9 C = (2d — l)/(2d) for N = 2. 
This becomes 5/6 for d = 3 (see Fig. []]). An interesting feature of the solution for 
the surface-anisotropy model is that the Curie temperature of the bilayer becomes 
independent of the lattice structure in the Ising limit r] s = 0. The result obtained 
above depends on the lattice dimensionality d only and, e.g., it is the same for the 
simple cubic model (d — 3) and the three-dimensional continuous-dimension model 
(d — 3.0) ||. In the Ising limit, the lattice structure comes into play for trilayers and 
thicker films, where the inhomogeneity of the gap parameter G n becomes essential. 

The thickness dependence of Curie-temperature shifts in films with small surface 
anisotropies are shown in Fig. ^ in the log scale. For N 1 they can be represented 
by the formula 



ttN c 3 k s 



k s = y/W{l/n. - 1), 



(2.25) 



with ^iuik 



Pz(l) — 1.51639 and the fitting parameter C3 ~ 1.36. This result, which 
will be derived analytically in the next section, is simpler than that for the model 
with the bulk anisotropy M, |J. The latter has the form (k c N <C 1) 

3 , 1 

lk 



r 1 

7 c,bv 



■ In ■ 

nN asK c N 



K c =y/2d(l/T}-l), 



(2.26) 



with 03 « 0.35, and iV under the logarithm makes the thickness dependence of 8 C shift 
substantially faster than 1/N (see Fig. ||). For n c N>l the bulk-anisotropy model 
shows a crossover to the finite-size-scaling regime described by Eq. (1.1) with A = 2. 
No such a crossover occurs for the model with surface anisotropy. 

In Fig. |^ the Curie temperatures of 100- and 200-layer ferromagnetic films with 
simple cubic structure are shown as function of the surface anisotropy. The film Curie 
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Fig. 3: Curie-temperature of ferromagnetic films with simple cubic structure vs surface 
anisotropy. 



temperature becomes greater than the bulk one for r/' s > 0.05. In this range it becomes 
independent of the film thickness, which is in accord with the surface character of the 
phase transition. Below the critical value of the surface anisotropy, the film Curie 
temperature falls below 6*k ulk . One can clearly see both the log dependence of the 9 C 
shift on the surface anisotropy and the 1/N dependence on the film thickness, as is 
given by Eq. (2.25). More careful analysis shows (see Sec. |]]) that the critical value of 
the surface anisotropy, which is defined from the condition T c (N,rj' s c ) = J^uik ^ en( j g 
to zero with the increase of the film thickness for d < 3. In three dimensions this 
dependence is logarithmic: rj' s C (N) ~ 1/lniV. This means that Eq. 2.25 is valid for 



sufficiently small anisotropy, r/' s <C rf B C (N), or, in other words, in the thickness range 
iVexp(— l/rj' s ) < 1. For whatever small value of rj' s , it will break down for very large 
N. Deviation of the numerically calculated points for rj' s = 10~ 2 in Fig. ^ downwards 
from the straight line corresponding to Eq. ( 2.26| ) is a manifestation of this incipient 
breakdown. 



3. Isotropic and weakly anisotropic films 

To get an idea about ordering in films with small surface anisotropics in d < 3 
dimensions, it is convenient to start with isotropic films. These films cannot order 
for any finite thickness N because they are systems of dimension d! < 2 and thus 
long-wavelength thermal fluctuations preclude ordering. On the other hand, it is 
physically clear that immediately below the bulk value of the Curie temperature the 
susceptibility of a thick fil m sho uld become extremely high. This means that the 
lowest eigenvalue \x\ in Eq. (2.13) becomes extremely close to zero. The contribution 
of this eigenvalue dominates in the the constraint relation (2.14), and this makes 
possible analytical calculation of [i\ . Since for the isotropic model there is no difference 
between longitudinal and transverse CFs, we will use here more compact notations 
E p and ip np [see Eq. ( [2.17 )]. First, Eq. ( 2.14 ) can be summed over all layers with the 
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use of the orthogonality of wave functions ^ n ip np ip np > = 5 pp i , which yields 



' 2d' „ / 2d' \ d'N 



^d' + Ep d ' \2d' + E p ) = ~W' ( ■ ' 

Next, one can subtract these equations for 6^ ulk and 9 from each other and separate 
the leading term with very small E\ (9) . This yields 

( 2d' \ d'N f 1 1 \ 

Pd ' {WTEW) )-^ N - — {9~9^J- (3 - 2) 

Since E p with p < 2 are expected to change not so strongly as E\ at the temperature 
interval 9^ ulk — 8, the quantity Eat can be expected to be subdominant in comparison 



to other parts of Eq. (3.2) 



For the simple cubic lattice, Pi is the Green function of the square lattice which 



is given by P 2 (X) ^ (1/tt) ln[8/(l - X)} for X = 1. Adopting this in Eq. fly) and 
exponentiating yields 



E 1 (8)^C N (9)exp 

where 



2irN (\ 1 



9 6>J? ulk 



(3.3) 



N 



C N (9) = E 1 {9^)X{^p^. (3.4) 

p=2 E PI?) 

Keeping high-lying eigenvalues with p ~ in the above formula is not justified, be- 
cause P2PO does not have its asymptotic form above in this case. On the other hand, 
the latter change negligibly for 9 close to d]? ulk and thus the corresponding numera- 



tors and denominators in Eq. (3.4) cancel each other. The low-lying eigenvalues also 



cannot change significantly in this temperature interval, thus the product in Eq. (3.4) 
should be of order unity. This leads to the order-of-magnitude estimation 

C N - Pi(^ ulk ) - C/N 2 , (3.5) 

which is sufficient for our purposes, since Cat will enter under the logarithm in the 
expression for the shift of the Curie temperature of the film. The second step in Eq. 
( |3.5| ) can be justified as follows. For thick films at the bulk criticality, G n is close 
to 1 in the main part of the film, excluding the regions near the surfaces. Thus for 
estimation of the eigenvalues one can set G n = 1 in the whole film, which amounts 
to the approximation i?i(6>J? ulk ) ~ £^ FA (6£ ulk ). Solution of the Schrodinger equation 



(2.16) with the potential V n — yields eigenvalues 



p 

[so that Ef FA (6£ ulk ) ~ 1/N 2 for N » 1] and eigenfunctions 

ipn P = C Np svc\(nq p ), C Np ~ 1/VN, (3.7) 

which describe a particle hopping in a rigid box. 

The picture described above is confirmed by numerical calculations, the results of 
which are shown in Fig. ^[ The latter were performed for the continuous-dimension 
model in d = 3.0. The dashed line in Fig. [| represents Eq. (3.3), where the transition 



E MFA ( 0bulk ) = 2(1 _ COS9p)) g = wp /(JV + l). (3.6) 



from the sc lattice to the d = 3.0 lattice is done by the replacement n 2, according 
to the general rule which can be found in Refs. [[| 0]. The constant C in Eq. ( |3.5| ) 
fits to 8 in d = 3.0 dimensions. 
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Fig. 4: Temperature dependence of three lowest eigenvalues E p for the isotropic film 
in d = 3.0 dimensions. Dashed lines represent Eq. (|3.3|) with tt =>■ 2 and cjv = 8/N 2 . 



One can ask how the variation of the gap parameter G n in the isotropic film 
below the bulk criticality looks like. The answer in the limit 9 <C 1 follows from the 
observation that all spins become strongly correlated and thus all a nn i become nearly 



the same for q = 0. Then from Eq. (2.2) immediately follows that b\ = &jv — \ and 
b n = 1 inside the film. This yields 



G r , 



2d/(2d- 1) 
1, 



n = 1,7V 
n ^ 1,7V 



(3.8) 



for 9 <C 1. The corresponding zero-temperature eigenvalues can be calculated analyt- 
ically and read 

E p (0) = 2(1 - cosq p ), q p = tt(p - 1)/7V. (3.9) 

These eigenvalues are all shifted downwards with respect to those of Eq. ([[!]), and the 
lowest eigenvalue is exactly zero, in accord with Eq. ( |3.3|) . The eigenfunction of this 
eigenvalue is constant throughout the film: ip n ,i = 1/yN. This is due to attraction 
of the particle to the potential wells at the boundaries of the box: Vi — Vn = — 1. 
Note that using Eq. (|j|) and Eq. (|J) in Eq. (Q yields C N = 0(1) at 6 < 1, in 
contrast to estimation (3.5) just below the bulk criticality. 



Calculation of the variation of the gap parameter G n in the film at 9 < 6^ ulk is an 



analytically intractable nonlinear problem, and the result of Eq. (3.3) does not help 
much. Linearization at 9 -C 1 shows that deviations of G n from the zero-temperature 
result of Eq. (3.5) are linear in temperature. A compact analytical solution can be 
only obtained for the trilayer. 

The deviation of the gap parameter from the bulk value, which is defined by 
Gi n = G n — G, is shown in Fig. || for the isotropic film in d — 3.0 dimensions at and 
slightly below the bulk criticality (in both cases G = 1). This deviation is proportional 
to the nonuniform part of the energy density || . At the bulk criticality, G n has the 
universal form 



G r . 



4 zK 

2dn 2 



d-3 



1 < n<7V/4, 



(3.10) 
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Fig. 5: Gap-parameter (energy-density) profile in the isotropic ferromagnetic film in 
d = 3.0 dimensions at and below the bulk criticality. 



for 2 < d < 4, as for the semi- infinite model jl0| [Tl], 0]. This yields the large-distance 
form 

V n = -{\- ii 2 )/n 2 , l«n<A/4, (3.11) 

for potential V n in Eq. ( |2.16| ). Note that in Refs. JlC], [□], the quantity V(z) = 
— (G n — 1) = 2dV n was used. For models with finite number of spin components, the 
energy-density profiles in critical films were calculated with the help of the e expansion 
||, |§. At T = 0.92T c bulk , the profile of G ln looks rather indefinite: in the middle 



of the film the tendency to the zero-temperature solution of Eq. (3.8) is already seen 



whereas closer to the boundaries G\ n still increases with lowering temperature. In 
the whole range of n, the relative deviation of G\ n from the bulk-criticality result is of 
order one. On the other hand, at such temperature the argument of the exponential in 
Eq. (3.3) is already —10, thus Ei is very small and futher lowering of the temperature 
leads to the instability of the numerical algorithm. Fortunately, the problem of finding 
the temperature variation of the gap-parameter profile in the film below the bulk Curie 
temperature becomes nonessential in the physically relevant three-dimensional case, 
because here the suppression of T c of the film is not so strong (see below). 

As we have see above, in isotropic films in d < 3 dimensions E\ is very small in 
a wide range of temperatures but turns to zero only at T = 0. If there is a s mal l 
anisotropy in the system, the basic equations for the transverse CF a nn ', Eqs. fl2.2j ) 
and (|]7), are slightly modified, and the variation of the gap parameter G n in the 
film slightly changes. These changes can be found perturbatively, although it is not 
easy to do analytically. When G n is inserted to the equation for the longitudinal CF 
<J Tm'i P er turbations of G n perturb, in turn, E\. Since Ei{9) goes almost horizontally, 
a small anisotropy is sufficient to cause Ei to cross the zero level at a transition 
temperature that is not small. 

The first step, finding the perturbed variation of the gap parameter G n , can be 
done qualitatively in the following way. If surface and bulk anisotropics, r)' s = 1 — ?? s 
and ?y' = 1 — 77, are small, one sets q = in the last term of Eq. ( 2.10| ), since this term 
creates a gap in the transverse CF a nn and it should be essential at small wave vectors. 
After that, defining Gn^ as the solution of the isotropic problem, one immediately 
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n/N 



Fig. 6: Numerically calculated squares of normalized eigenvectors ip np for the isotropic 
film (N ~ 500) in d = 3.0 dimensions at bulk criticality. Line ar behavior at small n is 



in accord with Eq. ( 3.15 ). The mean-field result of Eq. (3/7) is shown by the dotted 
line. 



finds that G n adjusts so that b n retains its isotropic value, i.e., 
d/( v G n ) + (1 - r,s/v)d'(S nA + 5 nN ) = d/G<®. 



(3.12) 



This defines the correction to G n due to anisotropy, which are positive. Now, pro- 
ceeding to the longitudinal CF, one can write for the eigenvalue problem of Eq. ( 2.1(\ ) 



V„ = K (0) +Vf J , where 

= -2d(l - V )/G<® - 2d'( v - Vs )(S nA + S nN ). 



r(i) 



(3.13) 



Numerical calculations show, however, that the surface part of this perturbation is 
somewhat oversimplified. It is not strictly localized in the boundary layer but re- 
distributed over some region, decaying in three dimensions slightly faster than 1/n 3 , 
presumably as l/(n 3 lnn). This feature in nonessential for the determination of the 
T c shift below; the difference of t he re sult with respect to those obtained with the 
use of the simplified form of Eq. ( 3.13| ) will be absorbed into analytically unknown 
numerical factors. 

The first-order corrections to E p due to the diagonal perturbation Vn have the 
form 



JV 



(3.14) 



as in the usual quantum-mechanical perturbation theory. In the temperature range 
of interest, slightly below the bulk criticality, the variation of the gap parameter G n 
does not strongly differ from that at bulk criticality. Thus estimation of ift^p m Eq. 



(3.14) can be done at T — . Here not too close to the surfaces and to the middle 
of the film one should consider the Schrodinger equation ( |2.16 ) with the potential 
V n given by Eq. (|1|). The standin g-wave solution of this equation for n > 1 in the 
semi-infinite geometry can be expressed through the Bessel functions and labeled with 
the wave vector taking continuous values from the interval (0, oo) |l0], |ll| §. In the 
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film this wave vector becomes quantized, and the normalized expression for tp np 
the region n < N/2 reads 



ipnp ~ {nqp/Ny^J^nqp), q p ~ irp/N. 



Thus for using J p {z) ~ z^ 1 , z <C 1, one obtains 



1+2 M , N 2(l+v)^ 



l/N 



2(1+M) 



(3.15) 



(3.16) 



For comparison, the MFA solution of Eq. (3.7) yields V'n ~ l/N 3 . Strong increase 
of the probability of finding the particle near the boundaries for N 3> 1 in our case, 
which is illustrated in Fig. |6[ is due to the long-range attractive potential V n . Note 
that at low temperatures, where G n approaches its limiting form given by Eq. (B.8), 
this effect becomes even stronger: (if, i — l/N. But there are no bound states near 
the surface in the is otrop ic model at any temperatures. 

Now, from Eq. ( 3.1 4| ) one obtains for the surface- and bulk-anisotropy models 



E 



(1, surface) 



-K 2 JN 



2(1+/*) 



E 



(l.bulk) 



(3.17) 



where in the bulk case G { n ] = 1 in the main part of the film and the normalization 
of eigenvectors ip np has been used. The Curie temperature of the film can be found 
from Eq. 
(O) and 



,(0) 



KTk) in the form E^' 
(\S-4)- Explicitly, one has 



Ei' 1 = 0, where E\"' is given for d — 3 by Eqs. 



r(0) 



{^Nf 



exp 



2ttN (1 



1 



(3.18) 



for t he su rface- and bulk-anisotropy models, respectively. This results in Eqs. ( [2.25 ) 
and ( |2.26 ), where the numbers C3 and 03 cannot be found analytically and should be 
fitted to the numerical solution. Remember that this analytical scheme for determi- 
nation of 9 C works if the argument of the exponential above is large. The method 
evidently breaks down for the bulk-anisotropy model, i f k c N > 1. Here the result for 
6 C crosses over to the finite-size-scaling solution of Eq. (1.1) with A = 2 and A ~ 1/k c 
H 0]. For the model with surface anisotropy, Eq. ( 3.18 ) also breaks down at suffi- 
ciently large N due to the surface phase transition, which will be considered in more 
detail in the next section. 



For very small anisotropy, the film orders at the temperature C <C 0, 



bulk 



where tpn 1 
yields 



l/N and in Eq. (3.3) Cn = 0(1) [see comment after Eq. (3 



-\N) 



ttN 



■In 



N 



-\N) 



-L ml 

ttN k c 



1, 

. This 
(3.19) 



for the surface- and bulk-anisotropy models, respectively, in accord with Eqs. ( |1.3| ) 
and (1.2). It should be stressed that the applicability conditions for the formulae 
above are difficult to fulfill for thick films, N 1 . For the l atter, the shift of the 
Curie temperature is typically small and Eqs. ( 2.25 ) and ( 2.26 ) are relevant. 

Let us consider now ferromagnetic films in dimensions lower than three. For the 
continuous-dimension model H the lattice Green function of layers, Pd', is given by 



Pd'{X) = -j-jT 



^dq 



1-XX' 



X' q - 1 - q 2 /(2d'). 



(3.20) 



For X close to 1 this yields 



C d >/K 2 d 7 d \ 



W dl + C d >K a d , 



d'-2 



d! < 2 
d! > 2, 



= y/2d'(l/X-l) < 1, (3.21) 
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where the Watson integral Wd> and the coefficient C'd> are given by 

K) 2 „ d' nd' 



{d'f 



c d > 



A d ' sin[(2 - d')n/2} ' 



(3.22) 



In three dimensions, the exact result is P3.0' = [1/(2A")] ln[(l + X)/(l — X)}. One 
should not mix up P3.0' (two continuous dimensions) with P2.0 (one discrete dimension 
and one continuous dimension), etc. 

For d < 3 the first term of Eq. (3.2) is of order 1/E± d ^ 2 and it dominates 
over Ejy. The latter is determined by other low- lying eigenvalues which are of order 
E p ~ (p/N) 2 . Thus Ejy ~ N 3 ~ d Cd'- This correction term will be retained in the 
formulae in order to provide correct limiting transition d — > 3. Using Eq. (3.21) and 



equating E\ to the anisotropic correction Ejp with the opposite sign, one obtains 



r 1 

7 c,bulk 



dC d , 
d'N 



(_#W)(3-d)/2 C d ' 



(3.23) 



Inserting here expressions for El 1 ' from Eq. (KVjh and using the value of p, from Eq. 
( 3.10| ), one arrives at the final results 



-1 

,bulk 



d' 



nd 



1 - {c d K s f- d N^ 2 l 2 



A d ' sin[(3 - d)n/2] (c d K s ) 3 - d N d - 2 +^ d ') 2 / 2 
for the surface-anisotropy model and 

d' Trd 1 - {a d n c Nf- d 



-1 

,bulk 



A d ' sin[(3 - d)w/2] 



nl- d N 



(3.24) 



(3.25) 



for the bulk-anisotropy model. Here q, Cd, and ad are numbers that should be fitted 
to the nume rical soluti on. O ne can check that for d — > 3 the formulae above go 
over to Eqs. ( [2.25 ) and ( |2.26| ) (the additional factors 2/ir in the latt er are due t o the 
differense between d — 3.0 and d — 3 models). Moreover, both Eqs. ( [3.24 ) and fl3.25| ) 
cross over to the single result Q~ x — 9~l uik ~ 1/N d ~ 2 [0 ^] in dimensions above 3, 
which is well-defined in the isotropic limit. 



4. Surface phase transition in films 

The surface of a semi-infinite magnetic system orders at a temperature above the bulk 
Curie temperature if there is a bound surface state of the Schrodinger equation ( [2.16 ) 



which lies below the continuum of the delocalized (bulk) states, i.e., E\ < Voo- In 
this case, with lowering temperature E\ reaches the zero value before all other (bulk) 
eigenvalues, and it dominates the longitudinal susceptibility \ zn — o^ n {<\ = 0)/$ 



[see Eq. (2.17)] in the boundary region, where the cigenfunction ip n i is localized. An 
example of the surface bound states is shown in Fig. [?] for the N = 8 film in d = 3.0 
dimensions at the bulk criticality. For the isotropic model, the potential wells near the 
surfaces are not strong enough to create a bound state. In contrast, for the extremely 
strong surface anisotropy the wells are deeper and there are bound states in each of 
the wells, which show a small tunnel splitting. Both models possess bulk states with 
E p > 0, which are not shown. The surface-anisotropy dependence of several lowest 
energy levels for thick films at the bulk criticality is shown in Fig. ||. One can see that 
the energy levels E p > nearly scale with 1/N 2 , which is characteristic for the bulk. 
Negative energies correspond to the states localized near the surfaces, here E\^ are 
nearly degenerate and practically independent of AT. 
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V„ = 2d(l/G H -1) 




1 2 3 4 5 6 7 „ 8 



Fig. 7: Numerically calculated effective potentials V n for the isotropic film (r)' a = 0) 
and that with the extreme surface anisotropy, rj s — 1, for A" = 8 in d = 3.0 dimensions 
at the bulk criticality. 



A natural question is how strong the surface anisotropy should be to create a po- 
tential well which can accomodate a bound state. The well known result in quantum 
mechanics (see, e.g., Ref. is that in one dimension a whatever small poten- 

tial well creates a bound state with the energy quadratic in the volume of the well: 

r m -I 2 

E = —4 V{x)dx for Ti = 1 and m = 1/2. If, however, the potential well is 

situated near a potential hump or a near an rigid wall, localization of the particle 
costs additional potential or kinetic energy, respectively, and it requires that the well 
strength exceeds some critical value. In this case the result is 

E = -A{P - P c ) c + const, (4.1) 

where P is appropriately determined strength of the well and for short-range potential 
wells C = 2. In the particular case of a rectangular well of depth Vq and width a, 
which is sided by a rigid wall, one has P = aq, q = y/Vo, P c — 7r/2, and A — ir 2 /(2a) 2 . 
If the potential V(x) has a long tail, the situation becomes more complicated, and 
the exponent £ deviates from 2, as we shall see b elow . 

Calculation of the critical strength P c in Eq. (4.1) for our problem ( [2.16 ) requires, 



in general, knowing the potential V n in the whole range of n including the surface 
region , n <~ 1. Whereas at the bulk criticality the asymptotic form of V n is given by Eq. 
( 3.11 ) , the values of V n for n ~ 1 can be only determined numerically || . Nevertheless, 



it can be shown that the isotropic semi-infinite model at the bulk criticality in d < 3 
dimensions is in the critical state, P — P c . A whatever small surface anisotropy r]' s 
makes the well deeper in the region n ~ 1 and it thus creates a surface bound state. 
As was argued above, this leads to the surface phase transition. This strong result 
follows from the form of the spin CF, which in the asymptotic region n,n' 3> 1 for 
q < 1 reads 0, 0, § 

0W(qJ - 2d0 < f £ (4.2) 



16 




with /i defined by Eq. ( 3.1 Cl| ) . Far from the boundary, qn, qn' ^> 1, this CF red uces to 



its bulk value, (T n n'{<\) = (d9 /q) exp(—q\n — n'\). In the region n,n' ~ 1 Eq. ( f4.2| ) is 
modified by nonsingular factors of order unity The spin CF above is proportional 
to the Green function which can be used to calculate perturbations of the solution 
of the problem (2.16), with E = q 2 , due to small perturbations of the potential V n . 



Such a perturbation theory, fails, however, in the bulk, since the bulk Green function 
above also diverges for q — *■ 0. A whatever small perturbation of V n changes the wave 
functions with E — > in a nonperturbative way, which leads to formation of bound 
states for attractive perturbations pjfl. To analyze the semi-infinite problem, one can 



use 

^^rraS" [1 + ° (4 z<<1 

^) = ^y[/-^)-W] (4-3) 

for the modified Bessel and Macdonald functions. One can see that for u < (i.e., 
d < 3) the Green function above diverges in the limit q — > (for d — 3 logarithmically), 
whereas for \x > (i.e., d > 3) it remains finite in this limit. Thus, in d > 3 dimensions 
there should be a critical value of the surface anisotropy, "q s ,c, above which there is a 
surface phase transition, whereas for d < 3 one has n s>c = 0. 

Different behavior for d > 3 and d < 3 observed above is entirely due to the 
different forms of V n for n *~ 1, whereas in the asymptotic region n» 1 the potential 
V n given by Eq. ( [3.11 ) is the same below and above three dimensions. If one goes away 



from d = 3 in both directions, the attractive tail of V n weakens, but for d < 3 the depth 
of the well increases in the surface region, n ~ 1, (see Fig. 1(6) of Ref. |(|), so that the 
well always remains in the critical state. In the limit d — » 2 the attractive tail of V n 
disappears, and the variation of the gap parameter G n approaches Eq. ([T^), which 
corresponds to V\ = Vn = — 1, V n = (n ^ 1, N). It can be checked directly that a 
whatever small further decrease of the boundary value of this potential leads, for the 
semi- infinite problem, N — oo, to the formation of a bound state. Determination of 
V n for n ~ 1 is an analytically intractable nonlinear problem. Nevertheless, Bray and 
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Moore |T^, [Tl| could obtain the spin CF of Eq. ( |4.2| ), which has different forms for 
d > 3 and d < 3 and contains the relevant information, without explicitly analyzing 
the region n ~ l! 

Now let us analyze how the energy of the surface bound state depends on the 
strength of the potential well if the latter slightly exceeds its critical value. For sim- 
plification, we will consider, instead of Eq. ( |2.16 ), a continuous Schrodinger equation 
—ip" + V(x)ifj — Eip with the potential V(x) modelled as 



V(x) = 




x < 
< x < 
x > a 



(4.4) 



[cf Eq. ( 3.11 )1- If we choose a — tt/2, then for d = 2 the long tail of V(x) disappears 
and Vo = 1 becomes the critical depth of the potential well, as for the original discrete 
problem. In general, this method tells nothing about the critical value of the surface 
anisotropy, but allows the determination of the exponent £ in Eq. (4.1). The bound 
solution of the problem above, if it exists, has the form 



i/j(x) 



C\ sin(<7a;), 

C2\fkxK Vi {kx\ x > a 



0<x<a {q = ^E + V ) 



(4.5) 



Here for very small \E\ one can neglect E in q and use the small-argument form of 
K^{z), which follows from Eq. ([i~g|). Then the boundary conditions at x = a result 
in the equation determining k: 



qa cot qa 



IMI 



1 — {kar^) 2 ^^ 



r(i-M) 



r(i + M) 



1/(2H) 



(4.6) 



Setting k — determines the critical value of the well strength P c , say, its depth Vo. 
For P slightly above P c , Eq. (|4.6|) can be represented in the form 



B(P - P c ) 



ImI(^v) 2 ^i ^ 

1 — {kar^) 2 \^\ ' 



which yields 



E 



-k 2 



BjP - P c ) 
B(P-P C 



d-3 



(4.7) 



(4.8) 



for the energy of the bound state. One can see that the "classical" one-dimensional 
behavior with the quadratic dependence of \E\ on P — P c is only realized for d = 2 
and d — 4 where \n\ = 1/2 and long tail of V(x) in Eq. (4.4) disappears. For d = 3 



Eq. (4.£) regularizes to the expression 



E ^ - (—\ exp 



B(P - P c ) 



7 = 0.5772, 



(4.9) 



which resembles the well known result for the energy of the bound state in two di- 
mensions [Q. Indeed, in two dimensions the radial part ip(r) of the wave function 
\f(r, 4>) = r~ 1 / 2 -0(r) exp(±im0), m = 0,1,2,... for the problem without potential 
energy satisfies the one-dimensional Schrodinger equation with the effective potential 
energy written in Eq. (4.4) for x > a, with /i m (see, e.g., Rcf. [19 ). Now, return- 
ing to the original problem with the surface anisotropy, one can notice that the depth 
of the potential wells near the surfaces change linearly with r]' s , thus one can replace 

P c by v' s 



in Eqs. (4.8) and (4.9) P 



rj' s c , where rj' s c = for d < 3. 
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Fig. 9: Critical surface anisotropy vs film thickness (50 < N < 500 for d — 3, 3.0, 3.5 
and 50 < N < 170 for d = 2.5) in different dimensions. The data correspond to the 
film with surface anisotropy on one of the two surfaces. The straight dashed lines are 
fits to the numerical data. 



The temperature of the surface phase transition, 6 Cl can now be determined using 
the results above. At 9 C , which is slightly above 6^ ulk , the energy of the bound state 
equals zero, but the bulk level of the potential, Voo, slightly exceeds zero. The value 



of 9 C can be found equating \E\, which is given by Eq. (4.S), to V"^: 



\E\ = Vat 



1-G 



nbulk\2z> b 



(4.10) 



where v\, = l/(d— 2) is the critical index for the bulk correlation length for the D = oo 
model. This yields 



nbulk 



~ V's.c) 



1 !/* 



ImI + b(v' s - V's.c) 



<l> = 



M-3| 

d- 2 



(4.11) 



The critical index <j> was calculated in Ref. j2(| for the model with arbitrary number 
of spin components n in the second order in e = 4 — d. In the limit n — > oo the result 
of Ref. ||20| becomes $ = l/2-e/4-£ 2 /8 + 0(e 3 ), which is in accord with Eq. ( |4.1l| ) . 



Note, however, that the e expansion fails below three dimensions for the model with 
infinite number of spin components, which is considered here. 

In films, surface bound states cannot be rigorously separated from the bulk ones. 
If these bound states are very shallow, which is the case near the special transition 
^surface _ gbuik^ ^ e i oca ii za tion length of the bound states is very large and it 
exceeds the thickness of the film. Because if this finite-size effect, the critical value of 
the parameter which controls the surface phase transition (here the surface anisotropy) 
cannot be determined unambigouosly. A natural choice is to define r] St c(N) from the 
condition 9 c (N,ri' s c ) — 6*J? ulk . The value of i] S7 c(N) can be found as the point of the 
intersection of the lowest eigenvalue E\ (N, r^) with the zero level at the bulk criticality 
(see Fig. ||). For the model with symmetric surfaces, however, the second eigenvalue, 
i?2, also goes down, crosses the zero level at somewhat larger value of -q' s and then 
very fast becomes almost degenerate with E\ . The latter situation corresponds to the 
two bound states well localized on both surfaces, with a small tunnel splitting. Thus, 
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crossing of E2 (N, r)' s ) with the zero level, as well as the degeneracy of E\ with E2 , could 
also be used as a criterion for the formation of bound states and thus for the special 
transition. One more and probably better possibility is to consider the film with a 
surface anisotropy on only one of the two surfaces. Here there is no complication 
arising from the tunneling between the bound states across the film; only the lowest 
eigenvalue E\ goes over to the bound state, whereas E2 remains always positive. 

In films there is no singularity in the dependence of E\ on the surface anisotropy, 
this dependence is linear near E\ — 0. On the other hand, the singularity of Ei(r]' s ) 
studied above for the semi- infinite problem above mirrors in the dependence n' r (N) . 
This dependence can be obtained if one uses a potential of the type of Eq. ( [4.4|) for 
a film, sets k = and imposes the symmetry condition on the wave function in the 



middle of the film, x = L/2. This yields Eq. (|4^7|) with P => P C {L), P c => P c (oo) 
and k =>■ 2/L. In terms of the original variables, dropping numerical factors, one can 
write 

- <o(oo) ~ ™ N _ m => ^ (d=3). (4.12) 

This result, as well as the conjecture t]' s c (oo) = for d < 3 made at the beginning of 
this section, are confirmed by numerical calculations the results of which are shown 
in Fig. I 

Positive values of <q' a C (N), even for d < 3, reflect the general tendency of the film 
to order at a temperature below the bulk Curie temperature. The latter is the case 
considered in the preceding section, and now it is clear that the applicability criterion 
for Eqs. (fT2f) and ( pi| ) is r]' s < rf a>c (N). For d < 3, Eqs. ( 0$ and ( p!^ ) break 



down for whatever small surface anisotropy rj' s , if the film thickness N is large enough. 
One can see that in three dimensions T]' s c {N) decreases logarithmically slowly, thus 



Eq. (2.25) works in a wide range N < exp{X/r]' s ) for small surface anisotropies. The 



applicability range of Eq. (3.24) shrinks fast with the descease of the spatial dimension 
d. 



5. Discussion 

In this paper it has been shown that the finite-size-scaling formula for the T c shift in 



magnetic films, Eq. (1.1), which seems to be the only theoretical tool for interpretation 
of experiments (see, e.g., Ref. [Q), describes in fact only one of several regimes. For 
the model with bulk anisotropy, the situation depends on ratio of the film thickness 
N and the transverse correlation length £ c j_ , which is usually ignored as a noncritical 
variable. For Nk c > 1, where k c = l/£ c i. at criticality, a different regime described by 



Eq. (2.26) H Q| is realized instead of Eq. (fL~l|). For the model with surface anisotropy, 



which is present in many experimentally investigated films, Eq. (1.1) never appears. 



Instead, the T c shift follows Eq. (2.25) in three dimensions for the surface anisotropy 



small enough. If surface anisotropy exceeds the critical value, rj' s > rf sc (N), Eq 



(2.25) breaks down and the film orders via the surface phase transition above the 
bulk Curie temperature (see Fig. |l|). A remarkable result is that rj' s C (N) goes to zero 
in the semi-infinite limit, N — > 00, for d < 3 (see Fig. ||). That is, a whatever small 
surface anisotropy leads to the surface phase transition in the bulk-isotropic semi- 
infinite model. This contrasts the isotropic model with enhanced surface exchange, 
which does not show any surface phase transition for d < 3. In three dimensions, 
f]' s C (N) ~ 1/lniV, thus Eq. ( 2.25[ ) is valid in a wide range of the film thicknesses: 



iV< exp(l/?^) for rj' s < 1. 

One can question whether all these effects, which have been demonstrated above 
for the D = 00 model, survive for the realistic classical Heisenberg model, D = 3. 
I expect that, in general, these effects should survive, because they are due to the 
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nearly Goldstone modes in a weakly anisotropic magnetic system, and these Gold- 
stone modes are inherent in all model with D > 2. On the other hand, the nonlinear 
coupling of fluctuations, which arises for finite number of spin components D, sup- 
presses fluctuations to some extent. This can be already seen from the fact that in 
the bulk 6 C = T C /T^ FA monotonically decreases with D and reaches its minimum in 
the spherical limit D = oo (strongest fluctuations). For the semi-infinite problem, 
the surface susceptibility xn a t the ordinary phase transition diverges for d < 3 (i.e., 
7°J d > for d < 3), if D = oo. For finite D the second-order e expansion (see, e.g., 
Ref. |Q and references therein, or, for a review, Ref. |2l]]) suggests that 7°J d remains 
positive at d — 3 (no divergence of the surface susceptibility) and probably changes 
sign at some critical dimension lower than 3. Thus, fluctuations are somewhat sup- 
pressed, and the situation is a bit closer to the mean- field one (d = 4), in comparison 
to the limit D = oo. This is an indication that in three dimensions a finite value of the 
surface anisotropy may be needed for the surface phase transition, in contrast to the 
zero value obtained in Sec. ^] Computing this critical value of the surface anisotropy 
with the help of MC simulations or other methods, as well as search for the regimes for 
the T c shift in films established above (or rather for their analogues for the Heisenberg 
model), seems to be an interesting problem. 
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